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The integrated Sachs- Wolfe (ISW) effect is an important imphcation for dark energy. In this paper, 
we have calculated the power spectrum of the ISW effect in the time varying vacuum cosmological 
model, where the model parameter /3 = 4.407 is obtained by the observational constraint of the 
growth rate. It's found that the source of the ISW effect is not only affected by the different 
evolutions of the Hubble function H{a) and the dimensionless matter density fim(a), but also by the 
different growth function D+ (a) , all of which are changed due to the presence of matter production 
term in the time varying vacuum model. However, the difference of the ISW effect in A(t)CDM 
' model and ACDM model is lessened to a certain extent due to the integration from the time of last 

^SJ \ scattering to the present. It's implied that the observations of the galaxies with high redshift are 

, ^ . required to distinguish the two models. 
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Q ! I. INTRODUCTION 

Since 1998, the type la supernova (SNe la) observations [l|, Q have showed that the expansion of our universe 
^] is speeding up rather than slowing down. During these years from that time, many additional observational re- 
I suits, including current Cosmic Microwave Background (CMB) anisotropy measurement from Wilkinson Microwave 
O ' Anisotropy Probe (WMAP) [1, |3], and the data of the Large Scale Structure (LSS) [l, [1], also strongly support this 
-^-^ : suggestion. In order to understand the mechanism of the accelerating expansion of the universe, the common recogni- 
^ tion has been received that there exists an exotic energy component with negative pressure, called dark energy (DE) , 
I— I , whose density accounts for two-thirds of the total energy density in the universe according to the current observations. 

It has been the most important thing to make an effort to probe into what the nature of DE is in the modern cosmol- 
' ' ogy. The simplest and naturalest candidate is the cosmological constant A (vacuum energy) 0, Q , with the equation 
I of state (EOS) w = ^ = —1. However, it's pity that it confronts with two difficulties: the fine-tuning problem and 
] the cosmic coincidence problem. An alternative offer is the dynamical DE to alleviate the difficulties above, such as 
' quintessence Q , phantom , quintom illj , generalized Chaplygin gas (GCG) model [H, Holographic DE Hi, and 
so on. 

At present, although in the presence of numerous dark energy models, the concordance model still fits best well with 
the current combined observations fl^ . So it's necessary to review the cosmic constant problems: The fine-tuning 
■ problem is the vacuum energy density predicted by quantum field theory (QFC) is lO^^Gey*, which is about 10^^^ 
orders of magnitude larger than the observed value 10~*^GeV', so there must be so-called "fine-tuning" while looking 
^ , for an balance between two numbers which exits so much discrepancy. The cosmic coincidence problem says: Since 
the energy densities of vacuum energy and dark matter turn on different evolutions during the expansion history of 
the universe, why are they nearly equal today and why it happens just now? To get this coincidence, it appears that 
. their ratio must be set to a specific, infinitesimal value in the very early universe. Therefore, in nature the second 
" " ' difficulty is also "fine-tuning" problem. In the final analysis, both of the two difficulties are related to the energy 
density of vacuum energy. On the basis of this consideration, a phenomenological and natural attempt at alleviating 
the above problems is allowing vacuum energy A to vary [8„il5,]. In A(t)CDM model, we keep unvaried for the EOS 
of the vacuum energy, equaling to —1, and have an time-evolving energy density. Although the weak point in this 
ideology is the unknown functional form of the A(t) parameter, we can deal with the A(i) parameter by different 
phenomenologically parameterized form, referring to Meanwhile, the underly explanation for A(t) parameter 

has also been explored in the framework of QFT using the renormalization group [l7| . 

The ISW effect [l^ arises from a time-dependent gravitational potential at late time and leads to CMB temperature 
anisotropies. From the physical picture point of view, the frequency of the photon released by CMB has been redshifted 
or blueshifted, which is caused by the decaying gravitational potential, while propagating along the line of light-like 
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from the last scattering surface to now. The ISW effect is a significant impUcation for dark energy. Since the 
gravitational potential was an constant during matter-dominated. At the late time, the gravitational potential is 
time-dependent when dark energy begin to dominate in the universe. So, the detection of ISW effect will provide 
great support and more information for dark energy. Now, the ISW effect is detected by the cross-correlation between 
CMB and LSS observations jT9l - l23j . In this way, we can separate CMB temperature anisotropics caused by ISW effect 
from the primary temperature anisotropics. 

In this paper, we consider the time varying vacuum model, where a parameterized form of A(t) is used. The 
densities of dark matter and vacuum energy don't conserve independently owing to the time varying vacuum term, 
which plays the role of the interaction between dark matter and dark energy. The matter production term leads to 
the different evolutions of the Hubble parameter H(a) and the dimensionless matter density fl^{a), and the different 
structure formation, comparing with the concordance model. Meanwhile, there being a relation between these three 
quantities and the source term of the ISW effect, we explore the impact from the change of these three quantities on 
ISW effect in the A(i)CDM model in theory. 

The paper is organized as follows. In next section, we briefly review of the cosmological evolution and structure 
growth in time varying vacuum model. In section III, we focus on the theory and detection of ISW effect and its 
power spectrum. The last section is the conclusion. 



II. REVIEW OF THE COSMOLOGICAL EVOLUTION AND STRUCTURE GROWTH IN TIME 

VARYING VACUUM MODEL 

Considering a spatially flat Friedmann- Robertson- Walker (FRW) cosmological model, which consists of two com- 
ponents: the dark matter and the dynamic DE. With the metric ds"^ = —dt^ + a?{t)[dr'^ + r'^{d6^ + sin^ Odcj)'^)], the 
Friedmann equation and energy-momentum conservation equation can be written as 

= ^^{Pm+ Pa), (1) 



'pi 

djpm + pa) 
dt 



= -3H{p„, + P„,+pa + Pa), (2) 



where, H = ^ is the Hubble parameter and Mpi = (SttG)^^ is the reduced Planck mass, p and P are density and 
pressure of a general piece of matter, and their subscripts m and A respectively correspond to dark matter and dark 
energy. 

In the present paper, we consider the time vacuum energy as dynamic dark energy with energy density pA — ^(O-^^p^ 
In order to simplify the formalism, we work in the framework of geometrical units with 8nG = c = 1. Thus, we can 
obtain the vacuum energy density p\ = A(t) and the corresponding equation of state Pa — -^PA — ^A(t). Then the 
equation (2) reads 

dpn 



dt 



3Fp„ = -A(i), (3) 



which shows that the decaying vacuum density A{t) plays the role of a source of matter production, and the densities 
of dark matter and vacuum energy don't conserve independently. Now, we consider a parameterized form using a 
power series of A(t) ^6|, i24|], this is : 

A{t) = mH + n2H\ (4) 
Combining the equations (1), (3) and (4), we can get the following Hubble function: 



Hit) - , (5) 



where we have defined /3 = 3 — n2. Then with the definition of Hubble function H{t) — ^, we proceed with the 
integration, deriving the evolution of the scale factor of the universe a(t): 

a(i) = ai(e^ (6) 

where ai is the constant of integration. Considering the two equations above, we obtain after some algebra that the 
Hubble function evolves with the scale factor of the universe as follows: 



ai 



(7) 
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where, the parameter ni and the integration constant ai can be reexpressed by taking the Hubble constant Hq and 
the current value of dimensionless dark matter density = 1) = ^mo into account with the help of equations (1), 

(4) and (7): 



ni = 



ai 



in 



mO 



niQ 



(8) 



Then making use of equation (8), we evaluate the dimensionless dark matter density flm{a) and the dimensionless 
quantity E{a): 



E{a) 



E{a) 1 - 



/3 



(9) 
(10) 



From equations (9) and (10), it's obvious that the Hubble parameter H{a) and the matter density ilm(a) scale 
differently, respectively comparing with their evolutions in the concordance model, due to the decaying vacuum 
density. In addition, we can get the comoving distance xi*^) from the solution of the Hubble function H{a) according 
to its definition of the integration ; 



Xia) 



da 



(11) 



In the following, we consider the structure growth in the framework of time varying vacuum model. The correspond- 
ing Newtonian equation governing the time evolution of dark matter perturbation could be generalized by taking a 
source of matter production into account [25|. It can be written as: 



D+ + {2H+ — )D+- 

Pm 



2 Prn V Prj 

5{x,a) 



D.=0, 



(12) 



where we have defined a growth function D+{a) = j^^^^^yy, an equivalent and convenient quantity describing the 

growth of dark matter density inhomogeneities. e is the source of matter production, and here it is given by e = — A(i). 
On account of the presence of e, it's seen that equation (12) is different from the well-known linear matter perturbation 
equation. 

There are two solutions to the differential equation (12), where one is a growing mode and the other is a decaying 
mode. Since almost all current models of the universe have a non-increasing Hubble rate, we are interested in the 
growing one, which is given by [igI [2^ : 



D+ia) 



a 2 



dx 



2/3 \ Silmo J Jo x'iE'^ix) 

Next, we constrain the A{t)CDM model parameter /3 according to the growth rate of clustering: 

a dD+(a) 



f{a) 



D+{a) da 



(13) 



(14) 



Using the currently available growth rate data, as shown in Table HI we can determine the best fit value of the 
model parameter /3 by minimizing 



x^(/?) = E 



fobs{Zi) - fth{Zi]l3) 



(15) 



where fobs is the observed growth rate, and cr/^^, is its corresponding la measurement uncertainty, fth is the theoretical 
value and can be obtained from equations (13) and (14). Here, we have taken a prior on 1)^0 = 0.27 [l|. The best fit 
value of model parameter /3 = 4.407 is obtained with the minimization iXmin/'i'^f — 2.83). In FIGlU we plot the 
observed values of the growth rate and its theoretical evolutions in ACDM model and A(i)CDM model with the prior 
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redshift z 


fobs 


Reference 


0.15 


0.51 ±0.11 


[26] 


0.35 


0.70 ±0.18 


[27] 


0.55 


0.75 ±0.18 


[28] 


0.77 


0.91 ±0.36 


[29] 


1.4 


0.90 ±0.24 


[30] 


3.0 


1.46 ±0.29 


[3^ 



TABLE I: The currently available observational data for the growth rate of clustering and its measurement uncertainty. 




FIG. 1: The observed values of the growth rate fobs with la uncertainty are plotted by the dots with error bars. The theoretical 
growth rate fth evolve with the redshift z for ACDM model (solid line) and for the A(f)CDM (dashed line). 
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FIG. 2: D+{a) vs. a for ACDM model (solid line) with n^o = 0.27 and A(f)CDM model (dashed line) with = 0.27 

and P = 4.407 on the left. The right panel shows dD+{a)/da scales with the scale factor a: the solid line is for ACDM 
model; the dashed line is for A(f)CDM model. The values of the cosmological parameter QmO = 0.27 and the Hubble constant 
Ho = 70.5km/s/Mpc are taken by [Ti ]. 



value of flmo = 0.27 and the best fit value of model parameter (3 — 4.407. From FIG.[1] it's found that the theoretical 
evolution of the growth rate in the ACDM model has a minor difference from the observational data at low redshifts. 

Then, the evolutions of growth function £'-|_(a) and its change rate dD^{a)/da as a function of the scale factor 
have been plotted in FIGlH From FIGlH we can see the growth function in A(t)CDM model is different from that in 
ACDM model for the presence of time varying vacuum energy. And it's found that is a little smaller than Da 

at a < 0.4 in the left panel, then there is a trend that the growth function in K{t)CDM model keeps larger than Da 
until now. In the right panel of FIG. [H it's found that when a < 0.16 and a > 0.65, the rate of change in the growth 
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function -D+(a) in A{t)CDM model is slower than that in ACDM model, but quicker than the change rate of Da in 
ACDM model with the scale factor ranges 0.16 < a < 0.65. And we can find that the increasing rate of the growth 
function Da (a) appears to be more and more small along all the scale factor. The right panel shows that there exists 
an upward trend in the increasing rate of the growth function D+(a) of A(t)CDM model till a = 0.3, then it gradually 
starts to decrease. In FIG[T]and FIG. [2 it's indicated that there is much discrepancy between the structure growth 
history in the two models. 



III. THEORY OF ISW EFFECT AND ANGULAR POWER SPECTRUM 



In this section, we give a brief review on ISW effect and its detection, and calculate its angular power spectrum. 
The ISW effect arises from a time-dependent gravitational potential at late time when dark energy dominates in the 
universe. The CMB temperature perturbation caused by ISW effect [sl] is expressed by 

2 / dT^ = -2 / dxa^Hia)^, (16) 



TcAiB J dr Jq da 

where r and x respectively the conformal time and comoving distance, which are related by dx = —dr. <I> 
denotes the gravitational potential, which is related to dark matter perturbation S according to Poisson equation on 
sub-horizon scales: 

V^$(,) = l^!^m, (17) 

Carrying out Laplace transform for the two sides of equation (17) and taking the dimensionless dark matter density 
rim(a) into consideration, we obtain: 

$(fc,a) = -^a^H'{an^{a)^-^. (18) 

Substituting the above equation into equation (16), we are able to write down the CMB temperature perturbation 
coming from ISW effect: 

AT ^rr2 . .daa)S{k^=l) 

= 3Ho / dxa H{a)— , (19) 

iCMB Jo da 

where the redefined function C(a) — a'^E'^{a)D,m{a)D^{a) is used. 

From equation (19), we can see that the source term of ISW effect depends on the following three quantities: the 
Hubble function H{a), the dimensionless matter density f2,„(a) and growth function D_|_(a), all of which are affected 
by the dark matter production term as we have mentioned in section II. 

From the observation point of view, ISW effect leads to these "new" temperature anisotropics, smaller than those 
produced at the last scatting surface. So there are many difficulties in directly detecting this effect apart from primary 
CMB temperature anisotropies. Now, this embarrassment is solved by searching for the cross-correlation between ISW 
temperature and LSS observations, which include the galaxies sample (ssl [33|. quasars sample jssj . and Luminous 
Red Galaxies (LRG)[36] from SDSS, the radio galaxy catalogue from the NRAO VLA Sky Survey (NVSS) [13, the 
infrared catalogue from two Micron All-Sky Survey (2MASS) [11], and the X-ray catalogue from the High Energy 
Astrophysical Observatory (HEAD) Here, we use the samples from SDSS and NVSS, which span different 

redshift ranges. The more detailed presentation on the characteristics of all the samples are given in Ref. [20l [2l|. 
The observed galaxies density contrast will be [13, HJ : 



Sg = j b{z)^Siz)dz = j f{z)5{z)dz, (20) 

where ^ is the selection function of the survey, h{z) is galaxy bias relating the visible matter distribution 5g to the 
underlying dark matter 5{z)^ which is either a constant or time-evolving with a parameterized form |40|. Here we 
adopt that it is scale independent for simplicity. 

We use the approximate function of the redshift distribution of the observed samples released by [20l [2ll l32l . mi, Si: 



7"! 

(21) 



(22) 
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where the parameters can be obtained according to the observation of the galaxy number density after we have selected 
a certain form of redshift distribution. For the SDSS galaxies sample with the redshift range 0.1 < z < 0.9, the best 
fittings of parameters in equation (19) are given by zq = 0.113, 7 = 1.197 and m = 3.457 with a median redshift of 
Zmed — 0.32 and a constant bias 6g = 1 in Ref. .20]. For the NVSS sample with a little wider redshift range, the 
best values of parameters in equation (20) are given by the effective bias beff — 1.98, = 0.79 and a — 1.18 in Ref. 
[2ll |. For the SDSS quasars sample with the redshift between 0.065 and 6.075, it's found that zq = 1.9, 7 = 2.2 and 
m = 2 and its median redshift Zmed ^ 1-8^^. With these best fit parameters, the normalized redshift probability 
distributions of the SDSS galaxies and quasars sample and NVSS sample are shown together in FIG|31 
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FIG. 3: the theoretical redshift distributions d N/d z of three samples with different redshift ranges. 
Introducing the weight function Wt{x} ^nd Wg{x) as follows [3^ : 

Wrix) = Sa^Hia)- 



da 
dN 



Wg{x) ^ H{a)h{a)—D+{a), 



(23) 
(24) 



the line of sight integral for ISW temperature perturbation I^T /Tcmb and the density contrast of observed galaxies 

H^S{k,a^ 1) 



Sg can be written as: 



AT r^" 

= / dxWTix)- ,2 

ICMB Jo k 

Sg= r" dxWgix)6{k,a^l), 



(25) 
(26) 



which allow the expressions for the ISW-auto spectrum Ctt{1), the ISW-cross spectrum Crgil) and the observed 
galaxies-auto spectrum Cgg{l) to be written in a compact notation, applying the Limber-projection 43] in the flat-sky 
approximation, for simplicity: 



Jo X^ 



k^ 



X 

Cggil) = I dx ' 



fc2 



Pss{k^l/X), 



(27) 
(28) 
(29) 



/o X 

where Pssik) is the present matter power spectrum. Here we take P{k) oc k"'T'^{k) and it can be normalized by 
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as = 0.8. T{k) is the transfer function, and we adopt its fitting form by |32l. |44|: 



T{q = k/ThMpc-^) = ^^^^^^^^'^'^^ [1 + 3.89g + {I6.2qf + (5.47g)3 + (6.71g)4] 



0.25 



(30) 



where the dimensionless quantity F = flmoh. 

lO-iOf 



10 



-11 



+ 



10" 



10 



-13 



10 



-14[ 



10 



15 20 



30 
1 



50 70 100 



FIG. 4: the ISW-auto power spectrum Ctt vs. multipole order I for ACDM model (solid line) and A(i)CDM model (dashed 
line). 
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FIG. 5: the ISW-cross power spectrum Cxg 



multipole order I for ACDi\f model: the short-dashed line is for the SDSS 



galaxies sample; the dashed line is for the NVSS sample; the long-dashed line is for the SDSS quasars sample; and A(t)CDM 
model: the black line is for the SDSS galaxies sample; the dotted line is for the NVSS sample; the dot-dashed line is for the 
SDSS quasars sample. 



As shown in FIG.|4l[5l[6l the power spectrums, including the ISW-auto power spectrum Ctt{1), the ISW-cross power 
spectrum Cxgil), and the observed galaxies-auto power spectrum Cgg{l), have been affected due to the presence of the 
time varying vacuum energy. The amplitudes of the ISW-auto power spectrums and ISW-cross power spectrums in 
A(i)CDM model are a few higher than those in the cosmological constant model. There is a minor difference between 
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10 15 20 30 50 70 100 

1 

FIG. 6: the observed galaxies-auto power spectrum Cgg vs. multipole order I for ACDM model: the short-dashed line is for 
the SDSS galaxies sample; the dashed line is for the NVSS sample; the long-dashed line is for the SDSS quasars sample; and 
A(f)CDM model: the black line is for the SDSS galaxies sample; the dotted line is for the NVSS sample; the dot-dashed line 
is for the SDSS quasars sample. 

the observed galaxies-auto power spectrum amplitudes in the A(t)CDMmodel and the ACDM model. The FIGH] 
and FIG. [5] show respectively the ISW-cross power spectrum CTg{l) and the observed galaxies-auto power spectrum 
Cgg{l) are related to the observed samples with different redshift ranges. In FIG. [5l it's found that the differences 
of the ISW-cross power spectrum Cxgil) in the two models become larger when the observed galaxies with wider 
redshift ranges are used. Since the total CMB temperature-auto power spectrum in ACDM model is still consistent 
with the observational results from WMAP, a few increase in the ISW-auto power spectrum Ctt{1) in A(t)CDM 
model may lead to the trouble between the total CMB temperature-auto power spectrum in A(t)CDM model and 
the observational results from WMAP. As we have mentioned, it's difficult to directly detect the amplitude of the 
ISW-auto power spectrum Ctt{1) apart from the total CMB temperature anisotropics. However, it's lucky that 
the ISW-cross power spectrum CrgH) provides an effective measurement for the ISW effect. In FIGl?! we show the 
observational data of the ISW-cross power spectrum Crg between ISW temperature and the radio galaxy catalogue 
from the NVSS and its theoretical evolutions in ACDM model and A(i)CDM model. It's found that the ISW-cross 
power spectrum Crg in A(t)CDM model is still consistent with the observational results within la errors. Although 
the amplitude of the ISW-auto power spectrum Ctt{1) in A(i)CDM model is a few higher than that in ACDM model, 
its contribution to the total CMB temperature power spectrum is still small and don't bring about the disagreement 
between WMAP observations and the total CMB temperature power spectrum. So, the WMAP measurements can 
not rule the A(i)CDM model out at present. In addition, FIGl5]shows the comparison of the observational data of the 
observed galaxies-auto power spectrum Cgg from the NVSS and the theoretical evolutions of Cgg in the two models. 

IV. CONCLUSION 

In summary, we have calculated the power spectrum of the ISW effect in the time varying vacuum energy model. 
It's shown that the source of the ISW effect is not only affected by the different evolutions of the Hubble function 
H{a) and the dimensionless matter density but also by the different growth function D+{a), all of which are 

changed due to the presence of matter production term in the time varying vacuum model. In FIG. [l]and FIG. [U 
it's seen that the difference between the growth function in A(t)CDM model and that in ACDM model is not minor. 
What's more, the left panel in FIG. [2] shows there is a crosspoint at a = 0.4 for the evolutions of the growth function 
in two models. The growth function in ACDM model is greater than that in A(t)CDM model before a = 0.4. When 
a > 0.4, the growth function in ACDM model is smaller than that in A(i)CDM model. However, since the ISW effect 
is calculated by integrating time-dependent gravitational potential from the time of last scattering to the present, the 
difference of the ISW effect caused by the difference of the growth function in the two models within the scale factor 
ranges [0.4, 1] can partly counteract that caused by the discrepancy in their growth functions at a < 0.4 though the 
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FIG. 7: Comparison of the observational data of the ISW-cross power spectrum Crg between ISW temperature and the radio 
galaxy catalogue from the NVSS [HI] and the theoretical evolutions of Cxg- the dots with error bars denote the observed data 
with its corresponding uncertainty; the solid line is for the ACDM model; the dashed line is for the A(f)CDM model. The 
well-measured temperature Tcmb is 2.725K by |14l. [45l|. 



0.00020 




FIG. 8: Comparison of the observational data of the observed galaxies-auto power spectrum Cgg from the NVSS [2ll | and the 
theoretical evolutions of Cgg: the dots with error bars denote the observed data with its corresponding uncertainty; the solid 
line is for the ACDM model; the dashed line is for the A(t)CDM model. 



integration. Therefore, the difference of the ISW effect in A(i)CDM model and ACDM model is lessened to a certain 
extent due to the integrated effect. As shown in FIG. El there is not much discrepancy between the availably observed 
quantities Cxg in the two models for the SDSS galaxies and the sample from NVSS. For the SDSS quasars with high 
redshift, the amplitudes of the ISW-cross power spectrum in A(t)CDM model and ACDM model obviously decrease. 
However, the relative discrepancy between ISW-cross power spectrums in A(t)CDM model and ACDM model for the 
SDSS quasars is larger than those for the SDSS galaxies and the sample from NVSS. Therefore, at present the ISW 
effect is weaker than the growth history in distinguishing ACDM model and A(i)CDM model. However, it's worth 
expecting the observational data of the ISW-cross power spectrum for the observed galaxies with high redshift in the 
future. 
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